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Abstract 

In this work, we calculate the form factors and the coupling constant of the strange-charmed 
vertex J/4’D*D S in the framework of the QCD sum rules by studying their three-point correlation 
functions. All the possible off-shell cases are considered, D s , D* and J/J>, resulting in three different 
form factors. These form factors are extrapolated to the pole of their respective off-shell mesons, 
giving the same coupling constant for the process. Our final result for the J/il)D*D s coupling 
constant is gj/^ DiDa = 4.30±o;3 2 GeV _1 . 
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I. INTRODUCTION 


In the past few decades, the Quantum Chromodynamics Sum Rules (QCDSR) community 
has been actively working in the calculation of strong coupling constants of a great variety 
of vertices. Particularly, our group has calculated a large number of charmed vertices as the 
D*Dtt [lj, D*Dp [2], J/ipD*D [3] among many others that were revisited in the review of 
ref. H|. Some recent works in QCDSR target vertices with the heavy beauty mesons ( B , 
B*. B s , B*, ...) as the B*BK |5j, B*B s i 7 ^ [ 6 ], Be* 1 B C Y f|[S ,/f ’ B,..I/v [ 8 ] and/or strange- 
charmed mesons ( D s , D*, D sl , D* s0 ,...) as the Jj^D s D s [21, TU] . J/ipD*D* [TU], J/^D* s0 D s i, 
<PD* s 0 D s1 HU, DlD a rP El and many others. 

Regarding the strange-charmed vertices, which are a subject of major interest nowadays, 
due to the detection of various new exotic states, it is noticeable that some of them have 
masses that lie just above the mass threshold. By exotic, we mean that the usual 

quark-model description as qq pairs does not hold for these new states. An example of 
such exotic state is the Y'(4140), first observed by the CDF collaboration [T2], which has an 
observed decay into the pair J/-00. Some interpretations for its internal structure are the 
tetraquark state ra, D S D* molecular state m, D*D* molecular state PHH] and hybrid 
state [T5] , Among them, we highlight the molecular D*D* hypothesis, so that the decay 
y(4140) —>■ can be understood with the help of an intermediate decay as U(4140) —>■ 

D*D* — > J/ipcj). In this decay, the vertex J/i(;D*D s is present, thus a more precise knowledge 
of its coupling constant and form factors may help in the understanding of the fundamental 
constitution of the Y'(4140) meson. Also, it may be useful in the studies of other exotic 
states as the A"(4350) [T 8 fj 2 T[ . F(4274) [22] and the F(4660) [23, 24], for example. 

In this work, we investigate the J/xj)D*D s vertex within the QCDSR formalism, following 
the development of our previous work for the J/^D S D S vertex [9]. This development consists 
in calculating the form factors from the three point correlation function in two ways: one 
with quarks and gluons degrees of freedom, the OPE side, and the other with hadronic 
degrees of freedom, the phenomenological side. After that, a double Borel transform is 
performed on each side and using the quark-hadron duality principle, we equate both OPE 
and phenomenological sides. If a window of stability is obtained, the form factor and coupling 
constant can be extracted. 

In the following sects., we present the QCDSR formalism and our results for this calcu- 
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lation. 


II. FORMALISM 

The starting point to perform the sum rule, is the calculation of the three point correlation 
function of the vertex. In this case, the J/p)D* s D s is a vector-vector-pseudoscalar (VVP) 
vertex. Each meson will be represented by an interpolating current, which contains the 
quantum numbers of the meson. 

We can calculate three correlation functions [3]: one with the vector meson J /ip off-shell, 
another with the vector meson D* off-shell and a third one with the pseudo-scalar meson 
D s off-shell: 

I < M \p :P ') = j d 4 xd 4 ye ip ' x e- iqy r^ M \x,y) (1) 

where q = p' — p is the transferred momentum, M is the off-shell meson (M = D s , J/ip, D*) 
and rW)(x, y) i s presented in eqs. @-0 below, 

r = < 0 ' I rof' (</) jf-'t ( 0 )} | 0 '), ( 2 ) 

T^\x,y) = {or\T{j?-(x)jy( V )j J /*'(0)}\0'), (3) 

(4) 

where T is the time ordered product and , j„ a and j$ s are the interpolating currents of 
the mesons J/ip , D* and D s respectively. 

According to the QCDSR, these correlation functions can be calculated in two different 
ways: using hadron degrees of freedom, called the phenomenological side, or using quark 
degrees of freedom, called the OPE side. In the following subsects., we calculate the corre¬ 
lations functions for the phenomenological and OPE sides, and the sum rule. 

A. The phenomenological side 

The effective Lagrangian, which represents the process in the J/ipD* s D s vertex, is [25ti26j: 

£j/iPd* s d s = -gj/yD* d s e^doppp (djD* s fi Dp + Dp) , (5) 

where e 0123 = +1 is the Levi-Civita totally antisymmetric tensor. 
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From this Lagrangian, we can obtain the vertices of the hadronic process necessary to the 
calculation of the phenomenological side of the QCDSR. For the J/-0, D* and D s off-shell 
cases they are, respectively: 

{D* s (p)J/^(q)\ D s(p')) = A)e*(p, A )q a p 1 £ a ^ 5 , (6) 

{J/i{>(p)D*(q)\D s (j/)) = ^i/;i* Ds (g 2 )^(p,A)e 5 (g,A)p a g 7 £ a/375 , (7) 

( J/il>(p)D a (q)\D*(p ')) = -igj^ D * Ds (q 2 )ep{p, A )e* s (p, A )p a p' 1 e aMS , (8) 

where (q 2 ) is the form factor of the J/ipD*D s vertex with meson M off-shell (M = 

J/il>,D s ,D* a ). 

In the calculation of the phenomenological side, we also make use of the following hadronic 
matrix elements: 


< 0 bf |P(g)> = (P(qM |0> = fp- m , p 7 , 

(9) 

( V (q)\jx\Q ) = fvm v e* x (q ), 

(10) 

(0| 3x\ v (<l)) = fvm v e x (q ), 

(11) 


where P is the pseudo-scalar meson (D s ), V is a vector meson (J/i i/j or D*), q is the four- 
momentum of the respective meson, m 2 PV is the squared mass of the meson, fpy are the 
decay constants, m qi and m q2 are the constituent quarks of meson P and e x is the polarization 
vector, with x = p(u) for (£>*). 

Finally, the correlation functions ([2])-(|4]) of the phenomenological side for each off-shell 
case read: 


= c 3j/*% D ,(.q 2 ) p'V^ . (12) 

^ (p 2 + m 2 D . ) (q 2 + m 2 ^) (p' 2 + m 2 Dg ) 

T phen(D* a ) = - C 9j^D t DSq 2 )P' X P^^Xa 

0 p 2 + m 2 J/ip )(q 2 + m 2 D ,) (p' 2 + m 2 Ds ) 

T phen(D s ) = C 9f/ S J Dt DM 2 )P' X P (T£ ^ hj . (14:) 

^ 0 p 2 + m2 j/^)(q 2 + m D a ) (p /2 + m D * s ) 

where h.r. stands for the contributions of higher resonances and continuum states of each 


-v2wA 


meson and C is defined as: 


C = 


f d s f d* fj/prn 2 Dg m Di mj 
(m c + m s ) 


(15) 


From eqs. (12)-(14), it is clear that this vertex has only one tensor structure to work 


within the formalism of the QCDSR. 
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B. The OPE side 


The OPE side is obtained from eqs. 0-0. with the interpolating currents written in 
terms of the quark fields: j J J v = cq^c, j^ s = cr) v s and j® 3 = ic^ 5 s. By construction, 
the OPE side is given by an expansion known as Wilson’s Operator Product Expansion, 
dominated by the perturbative term and followed by non-perturbative contributions: 

-pOPE(M) _ -ppert(M) _|_ -pnon-pert(M) (16) 


where pP® rt ( M ) j s the perturbative term and T™™ P ert ( M ) are the non-perturbative contribu¬ 
tions to the correlation function. When calculating the form factor, this expansion exhibits 
a rapid convergence and can be truncated after a few terms [1]. Considering the similarities 
between the J/il>D* s D s and both J/i)sD s D s and J/^>D*D vertices, we expect, from the for¬ 
mer, a similar behavior regarding the OPE series as seen in the two latter m- Therefore, it 
should be adequate to consider non-perturbative contributions up to the mixed quark-gluon 
condensate: 


-pnon-peri _ Wgg) , y*m q (qq) , ~p(p 2 G 2 ) , p( qgoGq) , -p7n q (qgaGq) 

ms ms I - ms ”T” ms ”T” ms ms 




p,v 


pis 


pis 


pis 


pis 


(17) 


Calculating eq. (16) using the non-perturbative contributions of eq. © corresponds to 
calculating the diagrams of fig. [lj As a consequence of the use of the double Borel transform, 


only the J/ifs off-shell case has contributions from all the non-perturbative terms of eq. (17). 
All the non-perturbative terms, except the gluon condensates (fig. 00. are suppressed 
for both D s and D* off-shell cases. This has been taken into account in fig. [lj wherefore 
these contributions were omitted from the beginning. The D s off-shell case is also omitted 
in fig. [lj as it is just an interchange between the mesons D* and D s of the D* off-shell case. 


Using dispersion relations, the perturbative term (fig. la) for a given meson M off-shell 
can be written in the following form: 


rr;‘ |M) (p.p') = 


POO vert(M) 
Pgv 


S, U, t ) 


47T 2 


(s — p 2 )(u — p' 2 ) 


dsdu , 


(18) 


where the spectral density pg e J t( ~ AI \s,u,t) can be obtained from the Cutkosky’s rules. The 
quantities s = p 2 , u = p' 2 and t = q 2 are the Mandelstam variables. 


In general, eq. (18) is the main contributing term of the OPE series in a QCDSR calcu¬ 


lation. For a VVP vertex, the spectral density can be parametrized as: 

= ^f (M) (s,U,l)/pV Act, 


(19) 
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where X = (u + s — t) 2 — Aus and F is an invariant amplitude. For the J/ijj and D s off 
shell, the invariant amplitude can be written as: 


— (m s — m c )(;A + B)—m s , 
p( D s) = _p( D t) — mc _ rn s )B — m c , 


( 20 ) 

( 21 ) 


where 


A = 


k 0 p' 0 |£;|cos# 
\p'\ y/s 


B = 


cos 6 = 


| A: | cos 6* 
\P'\ 

2p' 0 k 0 — u + (—l)h)(77 i 2 s — ml) 


2\p'\\k\ 


\ k \ = V k o ~ em2 s ~ (! - e )' m c , 


ko = 


s + e(m 2 s — ml) 


. s + u — t 
Po = 


2V~s ’ 

and e = 0(1) for Di*\j/ij>) off-shell. 


1^1 = 


a/A 
2y/~S ’ 


The first non-perturbative contribution to the correlation function (16) is the quark con¬ 
densate showed in fig. [Tb| 

-(ss)p ,X p a £^Xa 


r(ss)(J/4>) _ _ 

^ (p 2 - m 2 c ){p' 2 - m 2 c ) ' 


( 22 ) 


Figure lc is the first order mass term for the strange quark condensate contribution to 


eq. (16), numerically less important and usually negligible in the three point QCDSR: 


m c m s (ss) ip’ 2 +p 2 - 2ml) ,a 

2 P P ^fllP\(T ' 


pm s (ss)(J/'0) _ 

2 (p 2 — ml) 2 (p' 2 — ml ) 


(23) 


Full expressions for the gluon condensates (( g 2 G 2 )) of fig. Id li and mixed quark-gluon 
condensates ((sga ■ Gs)) of fig. lj lo for the J/if) off-shell case can be found in Appendix |Aj 


C. The QCD Sum Rule 

Two more steps are necessary in order to calculate the QCDSR to obtain the form factors. 
First, we make the change of variables p 2 —> —P 2 , p' 2 —> —P' 2 and q 2 —> —Q 2 followed by 


a double Borel transform [271 i28| to both sides of the QCDSR in eqs. (|12|)-(J14]) and ( |16| ), 
which involves the transformation: P 2 —> M 2 and P 12 —y M' 2 , where M and M’ are the 


6 




























Borel masses. After that, we equate the phenomenological and OPE sides, invoking the 
quark-hadron duality from which the sum rule is obtained. 


The second step is to eliminate the h.r. terms from the phenomenological side in eqs. (12) 


(14). This is achieved by the introduction of the continuum threshold parameters So and Uq 


in the OPE side and by taking advantage of the quark-hadron duality. These parameters 
fulfill the following relations: m? < So < m! 2 and rn 2 < uq < m ' 2 , where rn l and m 0 are the 
masses of the incoming and outcoming mesons respectively and m! is the mass of the first 
excited state of these mesons. 

After performing these steps, we obtain the QCDSR expressions for the form factors: 


%,m = 


_ 3 p rm ijrU / -i’> e -i^ e -Tfn,l s ,l rl _ 'y^-perr 
JSinf Ju in f y/X L J 


(Q 2 +™ 2 W 


9j/^d*dSQ ) — 


3 rs 0 ru 0 j_ F{Dl)e --^ e --^ dsdu + BB 

47T JUinf v A 


r<9 2 c 2 > 


c -m% s /M* 


(Q 2 +m) L) 


9j/Jd*d s (Q ) — 


— 2^2 f S ° f u ° ~rdsdu + BB 

4tt 2 J Sinf J U in f VA 


p(ff 2 G 2 > 


c 

(Q 2 +m 2 ) 


—m 2 . , /M 2 —7 

e e 


^ 2 / M ' 2 


where BB [ ] stands for the double Borel transform and BB [F non pert ] = BB 
the Di* off-shell cases, as already mentioned. The constant C is defined in eq. (15). 


T (a 2 G 2 ) 


(24) 


(25) 


(26) 


for 


As in previous works, the definition for the coupling constant gj/^D*D a is given by m 


lim 




(27) 


In order to calculate this, it is necessary to extrapolate the results for the form factor to the 
region of Q 2 < 0. From eqs. (24)-(26) it is clear that we can obtain the coupling constant from 


three different form factors, one for each meson off-shell. However, the coupling constant 
must be the same regardless the form factor used for the extrapolation. This condition is 
used to minimize the uncertainties in the calculation of the coupling constant, which will 
become quite clear in the following section. 
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III. RESULTS AND DISCUSSION 


Equations (24)-(26) show the three different form factors of the J/ijjD* s D s vertex. The 
numerical calculation of these form factors give results that must be fitted to an analytical 
function of Q 2 . In order to minimize the uncertainties regarding the fitting procedure, it is 
required that these three form factors lead to the same coupling constant |29j . This condi¬ 
tion also helps to reduce the errors of finding the Borel masses and continuum thresholds 
necessary to the calculation. 

2 

Regarding the Borel masses, which satisfy the relation M' 2 = M 2 , they can assume 

any value within the Borel window of stability, in which the phenomenologial and OPE 
side should be compatible. The Borel window can be reduced by requiring that the pole 
contribution must be bigger than the continuum contribution by at least 50% and that the 
perturbative term contributes with more than 50% to the total correlation function. We 
have worked with the average values of the form factors within the Borel window. Hence, it 
is not necessary to use any specific Borel mass value, minimizing the uncertainties associated 
to the choice of this quantity [9j. The procedure is briefly sketched as follows: the average 
value of the form factor is calculated in the Borel window for each value of Q 2 used. The 
standard deviation is then used to automatize the analysis of the stability of the form factor 
regarding the Borel masses and continuum threshold parameters. This stability criterion 
then indicates the optimal values of the continuum threshold parameters and the Borel 
windows we have to work with. Thus, not only is a good stability ensured in the Borel 
window but also in the whole Q 2 interval, without the need to chose a specific value for the 
Borel mass. 

The continuum threshold parameters, So and uq, are defined as sq = (m; + A*) 2 and 
Uq = ( m 0 + A 0 ) 2 , where the quantities Aj and A 0 have been determined imposing the most 
stable Borel window, as mentioned before. In order to include the pole and to exclude the 
contributions from higher resonances and continuum states, the values for A j/^, A D * and 
Ae> 3 cannot be far from the experimental value (when available) of the distance between 
the pole and the first excited state pm M\- Our analysis has found that the best values 
are Ad s = 0.6 GeV and Ad* = A j/^ = 0.5 GeV, which leads to remarkably stable Borel 
windows for the three off-shell cases, as it can be seen in fig. [2] for J/-0 and D s off-shell. 
Figures for D* off-shell are omitted as they are very similar to the D s off-shell case. 




Also in fig. [2j we can check that the perturbative term is the leading one in the OPE series, 
followed by the quark condensate (for the J/ip off-shell case) and the gluon condensates. The 
summed contribution of the non-perturbative terms, m s (ss), (sgaGs ) and m s (sgaGs) , is 
small (« 2%), which means we could easily neglect these terms without changing our results 
significantly. Likewise, the gluon condensates have a very small contribution to the total 
correlation function in the D s and D* off-shell cases and they may also be neglected. 

Besides the values of the Borel masses and the continuum thresholds, we also need to 
know the values of decay constants, condensates and hadrons and quark masses. The hadron 
masses used are m Ds = 1.968 GeV, mp* = 2.112 GeV and mj/^ = 3.097 GeV [50] . The 
other parameters are presented in table [TTJ 

In table |lj we present the Q 2 and Borel windows found for each form factor g^jf^ D * Da (Q 2 ) 
(where M is the meson off-shell), together with its parametrization and the corresponding 
coupling constant gj/VW Its associated error a is calculated using the method that will 
be explained later in this work. 

TABLE I: Parametrization of the form factors and numerical results for the coupling constant of 
this work. 


Quantity 

J/ip off-shell 

D* off-shell 

D s off-shell 

Q 2 (GeV 2 ) 

[0.5, 2.0] 

[1.0, 3.0] 

[1.0, 4.0] 

M 2 (GeV 2 ) 

[7.0, 7.9] 

[5.3, 7.3] 

[4.9, 5.9] 

9j^d*d s (Q 2 ) 

A 

TS+Cp 

Ae~Q 2 l B 

Ae~® 2 l B 

A 

193.4 GeV 

2.003 GeV 1 

2.330 GeV 1 

B (GeV 2 ) 

54.46 

6.027 

6.122 

9j/J>d*d s ^ a (GeV Q 

4 o 1+0.41 

4 20 +0,22 
^• zu — 0.18 

4 qo+0.27 
a.oj_ o 23 


The Borel windows presented in table[I](M 2 row) satisfy the already mentioned conditions 
regarding the pole and continuum contributions, as shown in fig. [3} 

The form factor obtained for the J/\b off-shell case was well adjusted by a monopo¬ 
lar curve, while for the D* and D s off-shell cases, the form factors were well adjusted by 
exponential curves (fig. [i]). 

The coupling constants of each off-shell case are shown in table [TJ They present different 
values among them, however, when the uncertainties are taken into account, they are all 
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compatible, as suggested by fig. [4j where the error bars stand for the uncertainties (a) of 
the coupling constants also presented in table |T[ 

The uncertainties were calculated following the same procedures used in our previous 
works mum- They are estimated by studying the behavior of the coupling constants when 
each parameter is varied individually between its upper and lower limits. Experimental 
parameters have their own errors (shown in table 0. while for the QCDSR parameters, we 
varied the thresholds A s and A u in ±0.1 GeV and the Q 2 window in ±20%. For the error 
due to the Borel mass (M 2 ), we used the standard deviation of the form factor within the 
Borel window. After that, we calculated the average of all these contributions and their 
respective standard deviations to obtain the coupling constant and its uncertainty for a 
given off-shell case. In table |TTJ we can see how the variance of each parameter affects the 
final value of the coupling constant. 

TABLE II: Percentage deviation of the coupling constants related to each parameter. 

Parameter_ A gj/$S;o,(%) A gj/jg;g,( % ) A 9j}$ DtDa (%) 


f Da = 257.5 ± 6.1(MeV)|&a 

1.94 

1.93 

1.94 

f D * = 301 ± 13(MeV) a 

3.53 

3.53 

3.53 

Jj/y = 416±6(MeV)[30] 

1.17 

1.18 

1.17 

m c = 1.27t°o° 0 l (GeV) [30] 

19.5 

12.9 

16.6 

rn s = 10lt 2 ? (MeV)[30] 

2.56 

0.96 

2.25 

M 2 (GeV 2 ) 6 

1.42 

4.46 

0.12 

± 0.1(GeV),A u ± O.l(GeV) 

8.93 

2.49 

2.08 

Q 2 ± 20% (GeV 2 ) 6 

4.94 

1.10 

0.98 

(ss) = -(290 ± 15) 3 (MeV 3 )[32] 

1.83 

- 

- 

(g 2 G 2 ) = 0.88 ± 0.16(GeV 4 )[33] 

1.45 

0.18 

0.73 

(sgcr ■ Gs ) = (0.8 ± 0.2)(ss)(GeV 5 )[33j 

0.49 

- 

- 

Fitting parameters ( A and B) 

14.9 

0.12 

0.06 


“Value obtained using /d*//d s = 1.17 from [351136] , 

6 The intervals for these quantities are those of table [l| 

According to table [HJ the variation of most of the parameters (the decay constants, 
condensates and the Borel mass for example) has a small impact on the value of the coupling 


10 







constant. On the other hand, the uncertainty in the quark charm mass is the one with the 
biggest propagation over to the final value of the coupling constant. The fact that the 
sensitivity related to the Borel mass is small was already expected, as we have very stable 
Borel windows (fig. [ 2 ]) . 


IV. CONCLUSION 


In this work we have calculated the gj/^D*D 3 coupling constant by three different QCD 
sum rules: one with the D s meson off-shell, another with the D* meson off-shell and a third 
one with the J/ij) meson off-shell. This procedure allowed us to reduce the uncertainties 
related to the method, leading to compatible coupling constants, as seen in fig. [4} 

Taking the mean value between the numbers presented in table |TJ we obtain the following 
final result for gj^ D * a D a '- 

gj/ipD* s D a = 4.30^q; 37 GeV -1 . (28) 


This coupling constant was obtained from sum rules that respect the criteria of the 
pole dominance over the continuum, the perturbative contribution of the OPE being the 
dominant one and the form factor stability regarding the Borel mass in the whole Borel 
window, as shown in figs. [2] and [3j 

Regarding the form factors, we observe some similarities with our previous work for the 
J/i\)D s D s vertex [SJ. Among such similarities, there is the form factor given by a monopolar 
parametrization when the heaviest meson {J/ip) is off-shell, while an exponential one is the 
case when one of the lightest mesons ( D s or D*) is off-shell. Furthermore, the OPE series 
presents, in both vertices, the same hierarchy for the contribution of each term, as well as 
comparable contributions among terms of the same dimension. 


We can also compare our result for the gj/^D* s D a coupling constant (eq. (28)) with the 


results of previous QCDSR works, presented in table III 


The comparison with ref. na is straightforward, its result differs from ours by ~ 30%, 
both being compatible within 2cr. Similar conclusions can be obtained from heavy quark 
effective theory (HQET) |39], from where the relation gj/^n*D s = gj/ii>D a D a l r ^D a stands, 


which allows us to compare the value of eq. (28) with the result of our previous work |9]. 
Using m 0j = 1.968 GeV, we have gj/ijjD a D a / r nD a = 3.04 GeV -1 , a different result (« 30%) 
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TABLE III: Values of coupling constants obtained using different methods. 


Coupling constant 

= 3.03 ± 0.62GeV _1 

-l 


Method and reference 

QCDSR Ill] 9j/ipD*D. 

QCDSR [ 3 j 9j/^d*d = 4.0 ± O.OGeV 

QCDSR i gj/y,D s D s = 5.98±g;«I 

VMD [371EH] gj/ijiDD = 7.44 


from the gj/^D*D s hi this work, but in agreement with the also sizable difference (~ 23%) 
found using the relation gj/^D*D = gj/ipDD /mo in [3j. A third comparison can be made 
invoking the SU(3) symmetry [35]: gj/^D*D s = gj/i/)D*D, from which we can compare the 


result of eq. (28) with the gj/^n*D of ref. [3]. From this symmetry, we can see that these 
two coupling constants are compatible with each other within lcr, with a difference between 
them of ~ 7%. This difference is comparable with the one found using the SU(3) relation 
gj/i!>D s D a = gj/ii)DD in ref. [9]. 

Finally, we can compare our result with the coupling constant gj/^DD from the vector 
meson dominance (VMD) [371 EH] using the SU(3) and HQET relations, which leads to 
gj/i>D*D a = 3.84 GeV -1 . This coupling is also compatible with the one of eq. (28) within 2cr. 
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Appendix A: Full expressions for (g 2 G 2 ) and ( qgoGq ) correlators - J/ij) off-shell case. 


Equation (Al) is the double Borcl transform of the gluon condensate contribution to the 
OPE. 


BB 


rV a ‘) 

(jr(P) 

967r 2 


r 


(n 2 G 2 ) 


x=d,e,f,g,h,i 

ai M , 2 t-t ( a l+ a lj£)[ m c- m s- a l m 2 c M ' 2 ] 

done M2 + Ml ' 2 “i m'' 2 - 1 p' x p a £^\v ^(oq), ( Al ) 

x=d,e,f,g,h,i 


where the index x = d to i corresponds to the diagrams of fig. Id li, with F^-Fp) given by 
eqs. (A2)-(AT) respectively. 

F(d) = m c {2a\m 2 c M 2 — 3aqM 2 + 2aqm c m s — 2aqm 2 + 3)/4 (A2) 

F( e ) = m c (a\tM 4 — 2 a 2 M 4 — 2 a\tM 2 + Aa\m c m s M 2 + 

6aqM 2 + apt + 4aqm 2 — 4aim c m s — 4)/(8(aqM 2 — 1)) (A3) 


F (f ) = (a\m c M 2 + m s - m c )/2 

F( g ) = 2m s (2a\m c m s M 2 — 3aqM 2 + 2aqm 2 — 2 aim c m s + 3)/(aqM 2 — l) 3 
= — (2aqm c M 2 + 3 m s — 2m c ) / (a\M 2 — 1) 

F(p = ( apm c M 2 + m s - m c )/(aq M 2 - 1) 


(A4) 

(AS) 
(A6) 

(AT) 


The double Borel transforms of the mixed quark-gluon condensate contributions are given 


by eqs. (|A8j) and (|A9j). 

B M iB M a [r< s "^ Gs >] = - S!in ^' S ‘. (3 [( Q 2 M' 2 M 2 + m 2 M 4 + 2 m 2 c M ,2 M 2 + m 2 c M' A ] - M 2 M 

J ~ Zj 1VI 1VJ 

2 2 

+3M 4 M /2 ) p /A p CT £^ A(T e - 2Pe~K'A (A8) 


B M 2B M ,2 [r m -< s ^ Gs >] = - ms ^ S9 ° G ^ (2 ml [ M 6 + 3M 4 M /2 + 3 M 2 M' A + M' 6 ] + 

VI 1VI 

2 2 

m c M 2 M' 2 [12M 4 + 6 M 2 M' 2 + M' 4 ]) p'^£^ Xa e~^e~^. (A9) 
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FIG. 2: OPE contributions for J/i^ off-shell (panel (a)) and D s off-shell (panel (b)). 
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FIG. 3: Pole and continuum contributions for the J/'ij) off-shell (panel (a)) and for D s off-shell 
(panel (b)), both at Q 2 = 1 GeV 2 . 
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FIG. 4: Form factors of the J/ipD*D s vertex. 
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